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Sizing  electrodes  is  an  important  step  during  Lithium-ion  battery  manufacturing  processes  since  poor 
cut  edge  affects  battery  performance  significantly  and  sometime  leads  to  fire  hazard.  Mechanical  cut¬ 
ting  could  result  in  a  poor  cut  quality  with  defects.  The  cutting  quality  can  be  improved  by  using 
a  laser,  due  to  its  high  energy  concentration,  fast  processing  time,  small  heat-affected  zone,  and 
high  precision.  The  cutting  quality  is  highly  influenced  by  operating  parameters  such  as  laser  power 
and  scanning  speed.  Thus,  we  studied  a  numerical  simulation  to  provide  a  guideline  for  achieving 
clear  edge  quality.  In  order  to  simulate  electrodes  laser  cutting  for  Lithium-Ion  batteries,  understand¬ 
ing  the  behavior  of  current  collectors  is  crucial.  This  study  focuses  on  current  collectors,  such  as 
pure  copper  and  aluminium.  Numerical  studies  utilized  a  3D  self-consistent  mathematical  model  for 
laser-material  interaction.  Observations  of  penetration  time,  depth,  and  threshold  during  laser  cut¬ 
ting  processes  of  current  collectors  are  described.  The  model  is  validated  experimentally  by  cutting 
current  collectors  and  single  side-coated  electrodes  with  a  single  mode  fiber  laser.  The  copper  laser  cut¬ 
ting  is  laser  intensity  and  interaction  time  dependent  process.  The  aluminium  laser  cutting  depends 
more  on  laser  intensity  than  the  interaction  time.  Numerical  and  experimental  results  show  good 
agreement. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  automotive  industry  is  strongly  motivated  to  satisfy 
consumer  demand  for  high  fuel  efficiency  and  environmental 
protection  agency  (EPA)  emission  regulation.  To  satisfy  these 
requirements,  there  are  efforts  to  accelerate  the  electrification  of 
future  powertrains.  The  market  for  hybrid  electric  vehicles  (HEVs), 
plug-in  hybrid  electric  vehicles  (PHEVs)  and  electric  vehicles  (EVs) 
is  rapidly  expanding,  especially,  in  the  US  and  Japan  [1  ].  To  develop 
more  efficient  hybrid  and  electric  vehicles,  an  understanding  of 
automotive  battery  technology  is  essential.  Among  the  many  bat¬ 
tery  technologies,  the  lithium-ion  battery  represents  the  most 
promising  technology  for  the  abovementioned  vehicles,  because 
of  their  high  energy-to-weight  ratio,  high  power  density,  a  lack  of 


Abbreviations:  L/V,  liquid-vapor  interface;  S/L,  solid-liquid  interface;  HEV, 
hybrid  electric  vehicle;  PHEV,  plug-in  hybrid  electric  vehicle;  EV,  electric  vehicle. 
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memory  effect,  and  lighter  weight  than  other  energy-equivalent 
secondary  batteries  [2,3].  The  lithium-ion  battery  consists  of  an 
anode,  cathode,  and  electrolyte.  The  anode  is  usually  graphite- 
coated  copper  and  the  cathode  is  lithium  metal  oxide-coated 
aluminum. 

One  of  the  important  factors  to  affect  the  performance  of  the 
lithium-ion  battery  is  the  manufacturing  process.  Electrodes  are 
produced  as  follows.  A  polymeric  binder  and  solvent  are  mixed 
with  active  electrode  powder  material  to  slurry.  A  thin  film  of  the 
slurry  is  spread  onto  both  sides  of  a  metal  foil.  Subsequently,  the 
coated  metal  foil  is  cured  in  a  drying  oven.  Then,  the  coated  mate¬ 
rials  are  cut  to  different  sizes  using  mechanical  tools.  Currently, 
die  cutting  and  rotary  knife  slitting  are  applied  to  size  the  elec¬ 
trodes  used  for  prismatic  and  cylindrical  cells,  respectively.  Both 
techniques  require  precise  and  relatively  expensive  tooling  that 
wears  over  time.  This  tool  wear  results  in  process  instability  and 
poor  cut  quality.  This  poor  cut  quality,  characterized  by  edge  bend¬ 
ing,  burrs  and  delamination  is  a  main  cause  of  short  circuits  [1] 
and  high  heat  generation  inside  the  battery  cell,  which  may  result 
in  the  catastrophic  failure  of  the  entire  module.  Laser  processing 
is  proven  in  the  industry  as  a  highly  efficient  and  reliable  manu¬ 
facturing  method  that  can  solve  these  problems  since  it  has  many 
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Nomenclature 

Physical  constant 

R  gas  constant,  8.314  x  1 015  J  K_1  mol-1 

o  Stefan-Boltzmann  constant,  5.67  x  1 0-8  W  m-2  K-4 

Symbols 

a  incident  angle,0 

A0  absorption  coefficient 

Cpi  liquid  constant-pressure  specific  heat,  J  kg-1  K-1 
CPs  solid  constant-pressure  specific  heat,  J  kg-1  K_1 
Cp  average  constant-pressure  mixture  specific  heat, 
J  kg-1  K_1 

Cpi  average  constant-pressure  specific  heat  of  liquid, 
J  kg-1  K-1 

Cps  average  constant-pressure  specific  heat  of  solid, 

J  kg-1  K_1 

8  delta  function 

ex,  ey ,  ez  unit  vectors  in  x,  y,  z  directions,  respectively 

s  emissivity 

fi  liquid  mass  fraction 

fs  solid  mass  fraction 

gi  liquid  volume  fraction 

gs  solid  volume  fraction 

F  force  function,  m  s-1 

Fevap  speed  function  of  the  L/V  interface  due  to  evapora¬ 
tion,  ms-1 

h  enthalpy,  J  kg-1 

hi  liquid  enthalpy,  J  kg-1 

hs  solid  enthalpy,  J  kg-1 

H£iC(0)  smoothing  function 

i  incoming  ray 

k  thermal  conductivity,  W  nrr 1  K_1 

ki  liquid  thermal  conductivity,  W  m-1  K-1 

ks  solid  thermal  conductivity,  W  nrr1  K_1 

k  curvature 

Lm  latent  heat  of  fusion,  J  kg-1 

Lv  latent  heat  of  vaporization,  J  kg-1 

ft  viscosity,  m2  s-1 

/JLyap  vapor  viscosity,  m2  s-1 

Hliq  liquid  viscosity,  m2  s_1 

m'eVap  net  mass  loss  due  to  evaporation  per  area, 
kgs-1  m-2 

n  normal  vector 

p  pressure,  Pa 

p  density,  kg  nrr3 

Pi  liquid  density,  kg  nrr3 

ps  solid  density,  kg  nr3 

q"  the  actual  net  energy  influx  at  the  liquid-vapor 

interface  per  area,  J  s-1  nr2 

the  spatial  laser  beam  distribution  after  multiple 
reflections  per  area,  J  s-1  nr2 

Qevap  energy  loss  due  to  evaporation  per  area,  J  s-1  nr2 
Rv  specific  gas  constant,  J  kg-1  K-1 

r  reflected  ray 

a  surface  tension,  N  m- 1 

t  time,  s 

T  temperature,  K 

Ts  temperature  at  surface,  K 

Tm  melting  temperature,  K 

Tb  boiling  temperature,  K 

Tv  temperature  at  vapor,  K 

u  fluid  velocity  vector,  m  s-1 

uv  velocity  of  vapor,  m  s-1 


x  spatial  vector 

0  level  set  function 


advantages,  such  as  contact-free  process,  high  energy  concentra¬ 
tion,  fast  processing  speed,  small  heat-affected  zone  (HAZ)  and 
flexibility  of  laser  power. 

The  contact-free  cutting  of  electrode  with  a  single  mode  fiber 
laser  for  the  lithium-ion  battery  has  been  investigated  to  deter¬ 
mine  the  achievable  cut  quality  [4,5].  A  comparison  between  a  CW 
mode  and  pulsed  mode  laser  was  studied  to  upscale  the  cutting 
speed  while  preserving  the  high  cutting  quality  [6].  It  is  important 
to  understand  the  physics  associated  with  laser  cutting  to  improve 
the  performance  of  the  battery.  To  the  best  of  the  authors’  knowl¬ 
edge,  it  has  not  been  studied  so  far.  To  fully  understand  the  behavior 
of  electrodes  during  the  laser  cutting,  investigation  of  the  behav¬ 
ior  of  the  current  collector  materials  is  an  essential  step.  Therefore, 
we  focus  on  the  current  collector  materials,  such  as  aluminum  and 
copper  in  this  study. 

Laser  cutting  involves  phase  changes,  such  as  solidification, 
melting,  and  evaporation.  A  Knudsen  layer  on  the  liquid-vapor 
(L/V)  interface,  where  the  continuum  hypothesis  fails,  is  formed 
due  to  a  high  evaporation  speed.  Moreover,  a  deep  penetration 
hole  is  observed  prior  to  accomplishing  the  cutting  as  shown  in 
Fig.  1,  which  requires  the  full  penetration  of  the  materials.  When 
a  deep  penetration  hole  exists,  laser  energy  can  accumulate  inside 
the  hole  so  that  the  absorbed  energy  increases  dramatically.  This 
increased  energy  absorption  might  in  turn  expedite  the  penetra¬ 
tion  hole  depth.  Due  to  these  physical  characteristics,  laser  cutting 
promises  cleaner  and  sharper  edge  compared  to  mechanical  cut¬ 
ting. 

Ki  et  al.  developed  a  mathematical  three-dimensional  self- 
consistent  laser-material  interaction  model  with  the  aid  of  the 
level-set  method  [7-12].  This  model  includes  fluid  flow,  heat 
transfer,  phase  changes,  multiple  internal  reflections,  free  surface 
evolution,  discontinuity  at  the  L/V  interface  and  surface  forces  on 
the  steel.  This  model  is  utilized  here  with  modifications  in  order  to 
investigate  the  physical  characteristics  during  laser  cutting  on  the 
current  collectors  of  the  lithium-ion  batteries. 

First,  this  paper  describes  the  mathematical  model  for  three- 
dimensional  self-consistent  laser-material  interaction  applied  to 
pure  copper  and  aluminum  foils.  Second,  the  simulation  results  are 
analyzed  and  discussed  taking  into  consideration  the  penetration 


Fig.  1.  A  deep  penetration  hole  during  laser  cutting  of  copper. 
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time  and  depth.  Based  on  the  numerical  analysis,  laser  parameter 
thresholds  are  provided  for  the  laser  cutting  of  the  current  col¬ 
lector  materials.  Third,  the  simulation  results  are  experimentally 
validated.  Finally,  the  conclusion  is  presented. 

2.  Mathematical  modeling 

2.1.  Governing  equations 


since  pure  materials  have  been  selected.  The  phase  enthalpies  are 
obtained  as 


hs  = 


h,= 


dT 

dT  +  Ljn  + 


/■ 

JTm 


CpidT 


(13) 

(14) 


The  three-dimensional  self-consistent  laser-material  interac¬ 
tion  model  is  summarized,  closely  following  Ki  et  al.  [  1 2  ]  with  some 
modifications  since  pure  aluminum  and  copper  foils  are  selected 
as  substrates.  This  study  assumes  an  incompressible,  laminar  and 
Newtonian  liquid  flow.  Vaporization  processes  are  treated  as  the 
volume  source.  Continuity,  momentum  transfer,  energy  transfer 
can  be  expressed  by  following  governing  equations  with  boundary 
conditions. 


where  Tm  is  the  melting  temperature,  Lm  is  the  latent  heat  of 
fusion,  Cpi  is  the  constant-pressure  specific  heat  of  liquid,  Cps  is 
the  constant-pressure  specific  heat  of  the  solid.  Furthermore,  the 
average  specific  heats  of  solid  and  liquid  (Cps  and  Cp/)  are  defined 
as 


dT 


(15) 


V  u=  emp8W) 

(1) 

+  V  •  (puu)  =  pV2u  “  “  ex  ■ 

(cm*:  —  vsr^)  S(4>) 

(2) 

+  V  •  (puv)  =  p,V2v  -  -  ey  ■  ^ 

crrnc  -  VST^  )  8(<t>) 

(3) 

9(g^}  +  V  ■  (paw)  =  mV2w  df  ez 

■  (oTlAT-  VsT^)  S(4>) 

(4) 

^LSe!I1  +  u  ■  V(pCp,T)  =  V  •  (kVT)  -  L 
+  q"(t,  x)5(0) 


d(pfi) 

dt 


+ 


d(pfsACpT) 

dt 


(5) 


where  u  is  the  liquid  velocity  vector,  x  is  the  spatial  vector,  p  is  the 
density,  /z  is  the  viscosity,  p  is  the  pressure,  k  is  the  thermal  con¬ 
ductivity,  Cpi  is  the  average-specific  heat  of  liquid,  t  is  the  time,  f 
is  the  liquid  mass  fraction  and  8(<f>)  is  the  delta  function  to  incor¬ 
porate  L/V  boundary  conditions.  Third  terms  on  the  right  hand  side 
of  Eqs.  (2)-(4)  are  boundary  conditions  for  the  L/V  interface.  This 
model  contains  the  capillary  and  thermo-capillary  forces  similar  to 
Ki  et  al.  [12].  The  second  and  third  terms  on  the  right  hand  side 
of  Eq.  (5),  which  are  the  boundary  conditions  for  energy  equations 
on  the  S/L  interface,  are  adopted  from  Bennon’s  continuum  model 
[13].  The  solid  mass  fraction  is  defined  as 


fs  =  1  ~fi  (6) 

Similarly,  the  volume  fraction  can  be  defined  for  liquid  and  solid 
(gi  and  gs).  The  relationship  between  the  mass  fraction  (//)  and  the 
volume  fraction  (fs)  can  be  described  as 

fs  =  (7) 

fl  =  ?f-  (8) 

P 

&+&  =  !  (9) 

With  these  relations,  the  density  (p),  thermal  conductivity  (/<),  and 
enthalpy  (h)  for  the  liquid  and  solid  mixture  are  defined  as 


P=gsPs+8iPi  (10) 


h  =fshs  +f,h, 


(12) 


J  1/77 

hs  =  CpsT 

hi  =  CpsTm  +  Lm  +  Cpi(T  —  Tm)  =  Cp/T  +  L 


(16) 

(17) 

(18) 


Cp  =  Cpl-fsCCpl-CpS)  (19) 

where  L  =  ( Cps  -  Cp/)Tm  +Lm.  By  using  the  defined  mixture  vari¬ 
ables,  the  boundary  conditions  for  S/L  interface  can  be  calculated 
[9].  The  energy  source  term,  which  is  a  fourth  term  on  the  right  hand 
side  of  Eq.  (5),  is  added  into  the  energy  conservation  equation  as  the 
boundary  condition  of  the  L/V  interface  in  Eq.  (5).  A  detailed  expla¬ 
nation  of  this  term  will  be  discussed  later  in  the  multiple  reflections 
section. 


2.2.  Level  set  equation  for  the  L/V  interface  tracking 


The  tracking  free  surface  plays  an  important  role  in  the  inves¬ 
tigation  of  the  interface  shape  and  the  process  physics.  The  level 
set  method  was  developed  by  Osher  and  Sethian  [14].  This  method 
has  increasingly  been  used  for  many  complex  problems  because 
of  following  reasons.  First,  it  is  easy  to  implement  and  conceptu¬ 
ally  straightforward.  Second,  surface  merging  and  separation  can 
be  automatically  handled.  Third,  geometric  qualities,  such  as  the 
surface  normal  and  the  curvature  on  the  surface,  are  easy  to  calcu¬ 
late. 

To  implement  the  level  set  method,  the  surface  of  interest,  which 
is  the  L/V  interface,  is  set  as  the  zero  iso-surface  or  zero  level  set. 
However,  the  equations  are  valid  for  both  a  zero  level  set,  as  well 
as  all  other  level  sets  [15]. 

0(x,  t)  =  ±d  (20) 


To  derive  PDE,  d  in  Eq.  (20)  is  set  to  zero  and  the  material  derivative 
is  taken  as 


D0 

Dt 


^  +  u  ■  V0  =  0 
dt 


(21) 


This  equation  can  be  re-written  as 


90 

dt 


+  u  •  n|V0|  =  0 


(22) 


Since  surface  normal  vectors  can  be  obtained  from  the  level  set 
function  as 


This  model  assigns  the  mass  liquid  fraction  1  and  0  when  the  tern-  n  _  ^ 
perature  is  above  and  below  melting  temperature,  respectively,  lv0l 


(23) 
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Fig.  2.  Smoothing  material  properties. 


Here,  u  n  is  the  speed  of  the  interface  in  the  normal  direction  and  F 
is  referred  to  as  the  speed  function  or  force  function.  Thus,  the  final 
form  of  the  level  set  function  is  a  hyperbolic  PDE. 

^+F|V0|=O  (24) 

The  speed  function  in  this  study  is  discussed  later  in  detail. 


2.3.  Smoothing  material  properties  near  the  L/V 


The  L/V  interface  has  significantly  different  properties  to  those 
of  the  S/L  interface.  Therefore,  the  material  properties  are  smoothed 
out  from  the  liquid  to  the  vapor  phase  using  discrete  functions 
proposed  as 


r  i 


HfilC(0)  =  < 


(p/  min(  A  z) 


if  4>  <  0, 
if  0  <  0  <  £, 
iif  0  >  £, 


(25) 


where  s  and  c  represent  the  width  of  the  interface  and  the  degree  of 
continuity,  respectively.  The  c  ranges  from  0  to  1 .  Thus,  the  material 
properties  can  be  expressed  by  this  function.  This  is  shown  in  Fig.  2. 
Viscosity,  for  example,  is  defined  as 


JLe,c  —  P'vap  +  (P'liq  ~  JLvap)Fls,C{(t)) 


(26) 


2.4.  Modeling  of  L/V  interface  phenomena 

In  the  laser-material  interaction  processes,  the  L/V  interface  in 
the  level  set  equation,  i.e.  0  =  0,  is  evolved  on  the  basis  of  fluid 
flow  and  vaporization.  The  fluid  flow  can  be  obtained  by  solving 
momentum  equation  proposed  in  Eq.  (3).  To  implement  vaporiza¬ 
tion  into  the  level  set  equation,  jump  conditions  with  back  pressure 
is  used  to  take  into  account  discontinuity  at  the  interface  cause  by 
the  Knudsen  layer  [16].  The  net  mass  flux  and  energy  flux  [7,16,17] 
can  be  obtained  as 


where 


\/2  RvTv 

F_  =  *Jmh[- 1  +  er/(m)]  +  exp(-m2) 

G  =  (2m2  +  1)[1  -  erf(m)]  — ^=m  exp(-m2) 

VVT 

_  2(2 m2  +  1  )^/Tv/Ts  -  2y/jtFn 
F-  +  ^TjfsG- 

Hence,  a  final  form  of  the  speed  function  is 


F  —  Fevap  +  u  •  n  (29) 

where  u  is  the  melt  pool  speed  and  n  is  the  normal  vector.  The 
vaporizing  mass  flux,  Fevap,  can  be  assumed  perpendicular  to  the 
L/V  interface. 


2.5.  Multiple  reflections  with  incident  angle 

The  energy  absorption  during  the  laser  cutting  process  is  greater 
than  the  energy  absorption  on  a  flat  material  surface.  The  portion  of 
total  energy  transferred  from  the  laser  beam  to  a  workpiece  can  be 
increased  dramatically  when  the  deep  penetration  hole  exists  due 
to  the  multiple  reflections  inside  the  walls  of  the  hole.  These  mul¬ 
tiple  reflections  can  be  identified  by  using  the  ray  tracing  method. 
If  the  laser  beam  is  irradiated  on  the  wall  of  the  deep  penetration 
hole,  then  the  reflected  rays  can  be  calculated  by 

r  =  i  +  2(-zn)n  (30) 

where  r  is  the  reflected  ray,  i  is  the  incoming  ray,  and  n  is  the  sur¬ 
face  unit  normal.  A  surface  unit  normal  vector  can  be  obtained  at 
each  point  with  the  aid  of  the  level-set  method.  In  this  method,  the 
laser  energy  is  absorbed  into  the  wall  at  each  reflection.  When  the 
rays  escape  the  computational  domain,  the  calculation  of  the  mul¬ 
tiple  reflections  is  terminated.  Furthermore,  incident  angle  is  taken 
into  account  in  determining  the  absorption  coefficient,  suggested 
by  Fabbro  and  Chouf  [18].  It  is  presented  as  following 

A(a)  =  Ao(cos  a)q  (31) 

where  a  is  the  incident  angle,  a  is  the  absorption  coefficient  on 
flat  surface,  and  A0  is  the  characteristic  of  the  material.  Again,  the 
incident  angle  can  be  calculated  based  on  the  information  from  the 
level-set  method.  The  values  of  the  material  characteristics,  such 
as  for  A0  and  q ,  vary  based  on  the  material  properties.  The  multiple 
reflections  contribute  a  large  portion  of  the  laser  energy  absorption 
in  the  deep  penetration  hole  during  laser  material  interaction  and 
it  increases  absorption  dramatically  compared  to  the  absorption 
on  the  flat  surface.  This  increased  energy  absorption  by  multiple 
reflections  causes  the  metal  to  behave  as  a  black  body.  After  con¬ 
sidering  energy  losses  due  to  evaporation  and  radiation,  the  actual 
net  energy  influx  at  the  interface  is 

=  Iv-Qevap-MiT4-^)  (32) 

where  is  the  spatial  laser  beam  distribution  after  multiple 
reflections  ct£(T4  -T^)  is  energy  loss  due  to  radiation,  a  and  e 
are  Stefan-Boltzmann  constant  and  emissivity,  respectively.  Note 
that  the  simultaneous  calculation  of  the  vapor  phase  and  the  liquid 
phase  allows  for  the  exclusion  of  convection  heat  loss  [12]. 


3.  Experiments 
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The  schematic  of  the  remote  laser  cutting  process  is  shown  in 
Fig.  3.  An  IPG  single  mode  cw  IR  fiber  laser  working  at  1070nm 
with  maximum  output  power  of  500  W  is  used  as  the  laser  source. 
The  laser  beam  is  fiber-delivered  using  a  1 0  p,m  core-fiber  diameter 
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Fig.  3.  Experimental  setup. 


and  is  collimated  to  a  beam  diameter  of  1 3  mm.  A  2D  galvo-scanner 
from  Scanlab  is  used  to  deflect  and  move  the  laser  beam  at  high 
dynamic  speeds  with  maximum  speeds  up  to  5  m  s-1 .  A  80  mm  F- 
theta  objective  mounted  on  the  scanner  keeps  the  scanning  image 
field  flat  within  a  35  mm  x  35  mm  square  field  size.  The  measured 
spot  size  at  the  focus  position  is  approximately  1 1  jxm  with  a  small 
Rayleigh  length  of  70  pan  and  M2  value  (86%)  of  1.3.  The  spot  is 
symmetrical  and  has  a  closely  Gaussian  energy  distribution.  Given 
the  small  Rayleigh  length,  it  is  essential  to  hold  the  material  flat 
on  the  fixture  during  the  entire  cutting  process.  A  vacuum  fixture 
is  used  to  hold  electrodes  and  prevent  any  movement  during  the 
cutting  process.  A  narrow  groove  (<1  mm)  is  machined  in  the  top 
plate  of  the  fixture  right  underneath  where  the  cutting  takes  place 
to  prevent  any  melt  attachments  that  come  out  of  the  cutting  kerf. 
The  fixture  is  mounted  horizontally  on  two  orthogonal  axes  of  a 
CNC  motion  system  while  the  scanner  is  mounted  on  its  third  ver¬ 
tical  axis.  This  precise  3-axis  positioning  is  required  for  the  setup 
to  align  the  fixture  with  the  scanner  and  to  position  the  scanner 
at  the  correct  vertical  distance  from  the  electrode  foil.  There  is  no 
process  gas  applied  during  cutting.  Cutting  is  performed  in  a  clean 
and  dust  free  environment  with  suitable  fume  exhaust  systems. 

Materials  used  for  laser  cutting  pure  copper  and  aluminum  foils 
with  a  thickness  of  10  pm  and  15[xm,  respectively.  In  addition, 
single  side-coated  electrodes  with  the  thickness  of  100  pan  are 
considered  to  validate  the  mathematical  model  in  a  more  realistic 
situation.  The  thickness  of  graphite  on  the  top  of  copper  is  90  pm. 
LiCo02  with  the  thickness  of  85  pm  is  coated  on  the  top  of  alu¬ 
minum.  The  electrodes  placed  on  the  fixture  as  the  active  electrode 
materials  face  upward  and  the  laser  beam  irradiates  on  the  surface 
of  active  electrode  materials. 

4.  Results  and  discussion 

4.1.  Simulation  preparations  and  solution  schemes 

The  numerical  domain  of  simulation  is  75  pm  x  30  pm  x  75  p,m 
with  the  non-uniform  mesh  size.  A  staggered  grid  is  used  to  obtain 
the  physically  proper  pressure  field  [19].  The  smallest  discrete 
mesh  distance  is  0.42857  pm.  Pure  copper  and  aluminum  are 


Table  1 

Material  properties  of  copper. 


Property 

Value 

Melting  temperature 

1357.77  (K) 

Normal  boiling  temperature 

2835.15  (I<) 

Critical  point  temperature 

8280  (I<) 

Liquid  density 

7920  ( kg  irr3) 

Solid  density 

8960  ( kg  irr3) 

Kinematic  viscosity 

3.50E-07  (m2  s1 )  [21] 

Surface  tension 

1.257- 0.0002  x  (7- 1356) 
(Nm4)  [22] 

Latent  heat  of  vaporization 

s^E+oeQkg-1) 

Latent  heat  of  fusion 

2.05E+05  (Jkg-1) 

Solid  thermal  conductivity 

317  (Wm-1  K-1)  [23] 

Liquid  thermal  conductivity 

157  (Wm-1  K-1)  [23] 

Liquid  constant-pressure  specific  heat 

571.6218  (Jkg-1  K-1) 

Solid  constant-pressure  specific  heat 

385  (Jkg-1  K-1)  [24] 

Liquid  thermal  diffusivity 

3.62E-05  (m2  s”1) 

Solid  thermal  diffusivity 

7.63E-05  (m2  s”1) 

Laser  absorptivity  for  flat  surface 

0.05 

Table  2 

Material  properties  of  aluminum. 

Property 

Value 

Melting  temperature 

933.47  (K) 

Normal  boiling  temperature 

2792  (K) 

Critical  point  temperature 

7963  (K) 

Liquid  density 

2333  (kgnrr3) 

Solid  density 

2700  (kg  m-3) 

Kinematic  viscosity 

4.43635E-07  (m2  s”1 )  [25] 

Surface  tension 

0.860-0.000115  x  (7-933.47) 

(N  m_1 )  [26] 

Latent  heat  of  vaporization 

1.09E+07  (Jkg-1) 

Latent  heat  of  fusion 

3.97E+05  0  kg-1) 

Solid  thermal  conductivity 

237  (Wm-1  K"1) 

Liquid  thermal  conductivity 

93.752  (Wm-1  K"1) 

Liquid  constant-pressure  specific  heat 

1255.2  (Jkg-1  K-1) 

Solid  constant-pressure  specific  heat 

896.9607116  (Jkg-1  K_1) 

Liquid  thermal  diffusivity 

3.20E-05  (m2  s”1) 

Solid  thermal  diffusivity 

9.79E-05  (m2  s-1) 

Laser  absorptivity  for  flat  surface 

0.07 

chosen  as  the  substrate  materials.  The  properties  of  these  materials 
are  shown  in  Tables  1  and  2. 

A  CW  single  mode  fiber  laser  with  a  focused  beam  diameter  of 
1 1  pm  is  used  and  it  is  moving  in  positive  x  direction.  A  Gaussian 
laser  beam  distribution  is  assumed  and  shown  in  Fig.  4.  The  grids 
for  multiple  reflection  calculations  are  much  finer  than  the  main 
grid  for  the  sake  of  a  good  resolution  of  solutions.  The  distance  of 
grid  for  the  laser  beam  is  0.1062  pun.  Each  mesh  has  different  laser 


Fig.  4.  Gaussian  laser  power  density  distribution. 
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Fig.  5.  Flow  diagram  of  the  solution  procedure  in  each  time  step. 


power  density  and  the  ray  tracing  technique  is  applied  to  the  laser 
beam  at  each  mesh  to  demonstrate  the  multiple  reflections.  The 
ranges  of  the  laser  power  and  scanning  speed  are  50-450  W  and 
1000  mm  s-1  to  5000  mm  s-1,  respectively.  Discretization  of  the 
governing  equations  is  accomplished  using  an  implicit  finite  differ¬ 
ence  method.  The  conjugated  Gradient  Stabilized  (CGSTAB)  method 
and  Semi-Implicit  Method  for  Pressure-Linked  Equation  Consistent 
(SIMPLEC)  are  adopted  to  solve  the  matrix  equation  and  coupled 
pressure-velocity  fields  on  the  staggered  grid  system,  respectively. 
Furthermore,  the  spatial  derivative  of  the  level  set  equation  is  dis¬ 
cretized  with  a  second-order  space  convex  scheme  [20].  Finally, 
the  value  of  c  for  all  of  the  material  properties  is  taken  as  1  and 
7  x  min(  Az)  is  chosen  for  e  in  this  study  to  smooth  the  material 
properties.  Az  is  the  smallest  mesh  size. 

A  flow  diagram  of  the  solution  procedure  in  one  time  step 
is  shown  in  Fig.  5  At  the  beginning  of  the  simulation,  initializa¬ 
tion  is  completed  to  tabulate  the  variables  on  the  Knudsen  layer 
and  to  assign  the  simulation  parameters.  After  the  initialization, 
multiple  reflections  are  executed.  The  temperature  and  velocities 
are  solved  sequentially  in  an  inner  loop  with  coupled  solutions 
until  their  computational  residuals  satisfy  certain  criteria.  Note 
that  a  maximum  iteration  is  set  to  prevent  unnecessary  compu¬ 
tations  due  to  slow  convergence,  even  though  the  residual  values 
may  not  be  satisfied.  After  obtaining  temperature,  velocity,  and 
pressure,  level  set  is  calculated  with  the  value  of  evaporation 
and  convection  that  affects  interface  deformation.  Finally,  physi¬ 
cal  properties  are  updated  with  a  given  liquid  mass  fraction.  This 
procedure  is  repeated  until  it  satisfies  the  designated  computation 
time. 

The  thicknesses  of  the  copper  and  aluminum  are  assumed  to 
be  10  pun  and  15  p,m,  respectively.  Variable  laser  parameters  are 
power  and  scanning  speed.  The  combinations  of  the  simulation 
conditions  are  tabulated  in  Table  3  for  copper  and  aluminum.  When 
the  depth  of  the  deep  penetration  hole  reaches  the  thicknesses  of 


Table  3 

Laser  processing  parameters  for  simulations  on  copper  and  aluminum. 


Sim# 

Copper 

Aluminum 

Power  (W) 

Speed  (mms-1) 

Power  (W) 

Speed  (mms  _1) 

1 

200 

5000 

100 

5000 

2 

250 

5000 

150 

5000 

3 

300 

5000 

200 

5000 

4 

350 

5000 

250 

5000 

5 

400 

5000 

300 

5000 

6 

450 

5000 

350 

5000 

7 

100 

3000 

50 

3000 

8 

150 

3000 

100 

3000 

9 

200 

3000 

150 

3000 

10 

250 

3000 

200 

3000 

11 

300 

3000 

250 

3000 

12 

50 

1000 

50 

1000 

13 

100 

1000 

100 

1000 

14 

150 

1000 

150 

1000 

15 

200 

1000 

200 

1000 

16 

250 

1000 

250 

1000 

Table  4 

Simulation  results  of  laser  cutting  penetration  time  and  depth  for  copper. 

Sim# 

Penetration 

Time  (|xs) 

Depth  (|jim) 

1 

1.8819 

0.957 

2 

6.7452 

4.787 

3 

0.9064 

10.000 

4 

0.5452 

10.000 

5 

0.4308 

10.000 

6 

0.3762 

10.000 

7 

0.0000 

0.000 

8 

3.2125 

0.531 

9 

6.7681 

1.808 

10 

4.3982 

10.000 

11 

0.8228 

10.000 

12 

0.0000 

0.000 

13 

0.0000 

0.000 

14 

0.7176 

0.106 

15 

6.8708 

2.233 

16 

3.5166 

10.000 

the  substrate  materials,  called  a  full  penetration  cutting,  penetra¬ 
tion  time  is  provided  following  simulation  termination.  When  there 
is  partial  cutting,  maximum  penetration  depth  during  the  simula¬ 
tion  is  obtained  first,  and  then  the  corresponding  penetration  time 
is  obtained.  Simulation  results  are  summarized  in  Tables  4  and  5. 

Table  5 

Simulation  results  of  laser  cutting  penetration  time  and  depth  for  aluminum. 

Sim# 

Penetration 

Time  (|xs) 

Depth  (|Jim) 

1 

0.4576 

1.382 

2 

0.3913 

15.000 

3 

0.2541 

15.000 

4 

0.2006 

15.000 

5 

0.1728 

15.000 

6 

0.1606 

15.000 

7 

3.5509 

1.382 

8 

0.4102 

1.382 

9 

0.4018 

15.000 

10 

0.2672 

15.000 

11 

0.2101 

15.000 

12 

3.1564 

1.382 

13 

8.5736 

1.808 

14 

0.3577 

15.000 

15 

0.2498 

15.000 

16 

0.2170 

15.000 
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Fig.  6.  Laser  cutting  penetration  time  for  copper. 
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Fig.  8.  Laser  cutting  penetration  depth  for  copper. 


4.2.  Thresholds  for  cutting 

The  laser  power  and  scanning  speed  thresholds  for  the  full 
penetration  cutting  are  obtained  as  shown  in  Tables  4  and  5.  For 
copper,  laser  power  between  250  W  and  300  W  is  the  threshold 
for  5000  mm  s-1  scanning  speed.  Laser  power  between  200  W  and 
250  W  is  the  threshold  for  both  3000  mm  s-1  and  lOOOmms-1 
scanning  speeds.  For  aluminum,  laser  power  between  100  W  and 
150W  is  the  threshold  for  5000,  3000,  and  lOOOmms-1. 

Laser-material  interaction  characteristics  can  be  derived  for 
both  materials  from  these  thresholds.  Interaction  time  and  laser 
intensity  are  closely  related  to  laser  speed  and  power.  Since  the 
thresholds  for  copper  vary  by  scanning  speed  and  laser  power, 
interaction  time  and  laser  intensity  are  both  effective  factors. 
However,  the  laser  power  thresholds  for  aluminum  are  the  same 
regardless  of  the  given  laser  scanning  speeds.  This  indicates  that 
the  aluminum  laser  cutting  is  more  affected  by  laser  inten¬ 
sity  than  the  interaction  time.  Therefore,  both  interaction  time 
and  laser  intensity  are  important  for  the  copper  laser  cutting. 


Moreover,  the  aluminum  laser  cutting  is  highly  influenced  by  the 
laser  intensity. 

4.3.  Penetration  time 

Penetration  times  of  full  penetration  cutting  for  copper  and 
aluminum  are  presented  in  Figs.  6  and  7,  respectively.  The  total 
penetration  times  during  full  penetration  cutting  for  copper  and 
aluminum  have  significant  differences.  The  penetration  times  for 
aluminum  are  always  less  than  1  p,s.  On  the  other  hand,  only  when 
laser  power  is  greater  than  300  W  is  penetration  time  less  than  1  p,s 
for  copper. 

For  copper,  the  penetration  time  decreases  exponentially  when 
laser  power  is  increased  from  250  W  to  450  W.  For  a  laser  power  of 
250  W,  no  full  penetration  cutting  is  obtained  with  a  5000  mm  s-1 
scanning  speed.  Moreover,  the  difference  in  penetration  time 
between  3000  mm  s-1  and  lOOOmms-1  scanning  speeds  is  larger 
than  other  laser  power  cases.  For  a  laser  power  of  300  W,  slight  dif¬ 
ferences  in  the  penetration  time  are  seen  among  the  5000,  3000, 
and  lOOOmms-1  scanning  speeds. 
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Fig.  7.  Laser  cutting  penetration  time  for  aluminum. 


Fig.  9.  Laser  cutting  penetration  depth  for  aluminum. 
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Fig.  10.  Threshold  and  kerf  width  based  on  simulation  results  and  experimental  results  of  laser  cutting  on  copper  -  dashed  line  is  the  kerf  width  of  the  simulation  (a)  250  W, 
5ms-1;  (b)  300 W,  5  ms-1;  (c)  200 W,  3  ms-1;  (d)  250 W,  3  ms-1;  (e)  200 W,  1  ms4 ;  (f)  250 W,  1  ms4. 


For  the  aluminum  case,  the  overall  trend  is  the  same  as  the 
copper  case.  However,  the  aluminum  case  shows  a  smoother 
decrease  of  penetration  time  with  increasing  laser  power.  While  the 
copper  case  exhibits  significant  differences  in  penetration  time 
with  the  laser  powers  of  250 W  and  300  W,  almost  no  dif¬ 
ferences  in  penetration  time  for  the  aluminum  case  are 
observed. 

These  results  show  again  that  the  copper  may  have  both  inter¬ 
action  time  and  laser  intensity  dependent  processes,  and  the 
aluminum  may  have  only  a  laser  intensity  dependent  process,  likely 
due  to  absorption  coefficients.  The  absorption  coefficients  for  a 
wavelength  of  1 .070  |xm  of  copper  and  aluminum  are  0.05  and  0.07, 
respectively.  The  higher  the  absorption  coefficient  of  the  mate¬ 
rial,  the  more  the  energy  that  can  be  deposited  on  the  material 
with  the  aid  of  multiple  reflections.  Therefore,  aluminum  is  pene¬ 
trated  more  rapidly  than  copper,  even  though  aluminum  foils  are 
thicker. 

4.4.  Penetration  depth 

Full  penetration  cutting  is  obtained  when  the  depth  of  the  pen¬ 
etration  hole  reaches  the  materials’  thicknesses,  which  are  10  [xm 
and  15  [xm  for  copper  and  aluminum,  respectively.  Once  full  pene¬ 
tration  occurs,  the  penetration  depths  are  automatically  assigned. 


However,  partial  penetration,  which  is  less  penetration  depth  than 
the  specimens’  thicknesses,  is  obtained  if  insufficient  energy  is 
provided.  Therefore,  a  relationship  between  the  laser  operating 
parameters  and  the  penetration  depths  could  be  observed.  These 
results  for  copper  and  aluminum  are  shown  in  Figs.  8  and  9, 
respectively.  Furthermore,  numerical  values  are  shown  in 
Tables  4  and  5. 

For  copper,  the  partial  penetration  depth  increases  expo¬ 
nentially  when  increasing  the  laser  power.  Values  of  the 
partial  penetration  vary  from  0.106  [xm  to  4.797  |xm,  which 
is  half  of  the  specimen  thickness.  On  the  other  hand,  for 
aluminum,  the  partial  penetration  and  full  penetration 
cutting  are  clearly  distinguished  between  100W  and  150W. 
The  maximum  partial  penetration  depth  is  1.808  |xm,  which  is  12% 
of  the  specimen  thickness.  From  these  results,  we  conclude  again 
that  the  aluminum  laser  cutting  is  laser  intensity  dependent. 

These  results  may  be  possible  due  to  the  materials’  proper¬ 
ties,  such  as  absorptivity  and  melting  points.  Since  aluminum 
has  a  high  absorptivity  as  well  as  low  melting  and  boiling  points 
compared  to  copper,  aluminum  absorbs  more  energy.  Once  alu¬ 
minum  absorbs  enough  energy,  the  material  is  evaporated  at  a 
relatively  fast  speed.  Therefore,  the  aluminum  laser  cutting  shows 
clear  differences  between  partial  penetration  and  full  penetration 
cutting. 
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Fig.  11.  Threshold  and  kerf  width  based  on  simulation  results  and  experimental  results  of  laser  cutting  on  aluminum  -  dashed  line  is  the  kerf  width  of  the  simulation  (a) 
100W,  5ms-1;  (b)  150W,  5 ms”1 ;  (c)  100 W,  3 ms”1 ;  (d)  150W,  3 ms”1 ;  (e)  100W,  1  ms4;(f)  150W,  1  ms4. 


4.5.  Experimental  validation 

A  series  of  linear  laser  cutting  tests  of  the  current  collector  mate¬ 
rials  and  single  side-coated  electrodes  are  performed  for  validation 
with  combinations  of  the  laser  parameters,  as  shown  in  Table  4. 
Here  the  laser  power  is  the  calibrated  power  output  on  the  work- 
piece  and  includes  all  the  beam  delivery  losses.  During  all  cutting 
tests,  the  focus  of  laser  beam  was  positioned  on  the  surface  of 
the  electrode  foil  material.  The  length  of  cut  was  fixed  at  25  mm 
and  only  a  single  pass  of  laser  beam  was  used  for  cutting.  The 
start  and  stop  regions  of  the  cut  length  were  not  considered  for 
analysis  to  avoid  any  effects  of  scanner  acceleration  and  deceler¬ 
ation  at  high  speeds.  To  ensure  accuracy  of  test  results,  cutting 
tests  were  repeated.  Among  these  experimental  results,  SEM 
analysis  is  done  for  the  cases  that  have  laser  power  and  laser 
scanning  speeds  near  to  full  penetration  cutting  thresholds.  The 
threshold,  which  is  highlighted  area,  and  top  view  of  the  copper 
laser  cuttings  are  shown  in  Fig.  10,  using  a  different  combi¬ 
nation  of  laser  parameters.  The  kerf  widths  obtained  from  the 
aluminum  laser  cutting  simulations.  The  kerf  widths  of  the 
experimental  and  simulation  results  show  a  good  agreement 
for  the  full  penetration  of  high-speed  cutting  cases  Fig.  10(a) 


shows  almost  no  cutting.  The  small  area  of  full  penetration 
is  barely  seen.  Reattachments  in  the  cut  surface  are  shown  in 
Fig.  10(c).  A  bubble-shape  recast  is  observed  with  inconsistent 
kerf  width  in  Fig.  10(e).  This  bubble-shape  recast  creates  the  bad 
cut  surface  of  the  electrodes  so  that  battery  performances  could 
be  deteriorated.  On  the  other  hand,  Fig.  10(b)  and  (d)  shows 
clear  cutting  and  consistent  kerf  width.  Even  though  250  W  and 
1000  mm  s-1  laser  beam  create  the  full  penetration,  bubble-shape 
recast  is  seen  in  Fig.  10(f). 

Fig.  1 1  shows  the  threshold,  which  is  highlighted  area,  and  top 
view  of  the  laser  cutting  of  aluminum.  Dashed  lines  depict  the  kerf 
widths  obtained  from  the  simulation.  The  kerf  widths  of  the  exper¬ 
imental  and  simulation  results  show  a  good  agreement  for  the  full 
penetration  cases.  Fig.  11(a)  shows  reattachment  at  the  cutting 
zone.  These  reattachments  bridge  the  two  separated  foils.  Depend¬ 
ing  on  the  size  of  this  reattachment,  some  of  them  dissociate  easily 
when  merely  touched.  Similar  phenomena  could  be  observed  in 
Fig.  11(c).  Thus,  we  could  treat  these  cases  as  the  partial  penetra¬ 
tion  cutting.  On  the  other  hand,  full  penetration  cutting  is  achieved 
in  Fig.  11(b),  (d),  and  (f),  which  are  150W  laser  power  cases,  as 
predicted  by  the  simulation  results.  However,  Fig.  11(e)  shows  full 
penetration  cutting. 
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Fig.  12.  Threshold  and  kerf  width  based  on  simulation  results  and  experimental  results  of  laser  cutting  on  anode  -  dashed  line  is  the  kerf  width  of  the  simulation  (a)  100  W, 
5ms-1;  (b)  150 W,  5 ms-1 ;  (c)  100 W,  3  ms”1 ;  (d)  150 W,  3  ms”1 ;  (e)  100 W,  1  ms4;(f)  150 W,  lms4. 


Therefore,  simulation  results  match  well  with  experimental 
studies  for  the  high  speed  laser  cutting,  such  as  5000  mm  s-1  and 
3000mms_1,  of  the  current  collector  materials.  However,  there 
is  a  discrepancy  for  the  low  speed  cutting,  such  as  lOOOmms-1. 
The  reason  for  this  might  be  that  the  mathematical  model  used  a 
constant  absorption  coefficient.  Therefore,  including  a  temperature 
dependent  absorption  coefficient  might  improve  the  mathematical 
model  for  low  speed  cutting. 

The  kerf  widths  and  thresholds  of  current  collector  materi¬ 
als  obtained  by  simulation  are  compared  to  the  experimentally 
obtained  results  of  current  collector  materials  in  the  presence  of 
the  active  electrode  materials  to  validate  the  proposed  mathemat¬ 
ical  model  in  a  more  realistic  situation.  The  threshold  of  copper 
laser  cutting  and  top  view  of  the  single  side-coated  anode  laser 
cuttings  are  shown  in  Fig.  12.  Only  graphite  is  sublimated  in 
Fig.  12(a).  Fig.  12(b)  shows  the  partial  cutting  of  copper.  Thus,  the 
laser  power  of 300  W  provides  not  enough  energy  to  obtain  full  pen¬ 
etration  cutting  with  the  5000  mm  s-1  scanning  speed.  With  the 
3000mms_1  scanning  speed,  simulation  results  match  well  with 
experimental  studies  as  shown  in  Fig.  12(c)  and  (d).  Although  the 
kerf  width  of  copper  in  Fig.  12(d)  shows  a  good  agreement,  wider 
kerf  widths  are  observed  for  graphite.  Both  Fig.  12(e)  and  (f)  show 


full  penetration  cutting  with  3.7  times  wider  kerf  widths  compared 
to  simulation  results. 

Fig.  13  shows  the  thresholds  of  aluminum  laser  cutting  and  top 
view  of  the  single  side-coated  cathode  laser  cuttings.  No  cutting 
is  observed  in  Fig.  13(a).  Regardless  of  laser  scanning  speeds,  the 
laser  power  of  150W  provides  full  penetration  cutting  as  seen  in 
Fig.  1 3(b),  (d),  and  (f).  This  full  penetration  is  predicted  by  the  simu¬ 
lation  results.  Furthermore,  the  kerf  widths  of  aluminum  obtained 
experimentally  and  numerically  are  in  good  agreement,  as  shown 
in  Fig.  13(b),  (d),  and  (f).  However,  Fig.  13(c)  and  (e)  shows  full 
penetration  cutting  and  these  are  different  to  results  estimated  by 
simulation. 

Simulation  and  experimental  results  of  copper  with  the  sin¬ 
gle  side-coated  graphite  are  in  good  agreement  in  the  case  of  the 
3000mms_1  scanning  speed.  In  addition,  simulation  results  of 
aluminum  with  single  side-coated  LiCo02  match  well  with  exper¬ 
imental  studies  in  the  case  of  5000  mm  s-1.  The  other  cases  show 
discrepancies.  These  discrepancies  can  be  caused  by  the  differ¬ 
ent  absorption  coefficients  between  current  collector  and  active 
electrode  materials.  Furthermore,  the  sublimation  characteristic 
of  graphite  might  affect  these  discrepancies.  Finally,  the  composi¬ 
tion  variation  of  current  collector  materials  between  the  interface 
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Fig.  13.  Threshold  and  kerf  width  based  on  simulation  results  and  experimental  results  of  laser  cutting  on  cathode  -  dashed  line  is  the  kerf  width  of  the  simulation  (a)  100  W, 
5  ms-1;  (b)  150  W,  5  ms-1 ;  (c)  100  W,  3  ms”1 ;  (d)  150  W,  3  ms”1 ;  (e)  100  W,  1  ms4;(f)  150  W,  lms4. 


of  current  collector  and  active  electrode  materials  may  result  in 
discrepancies  between  simulation  and  experimental  results. 

5.  Conclusions 

The  numerical  studies  of  the  mathematical  model  have  been 
described  for  cutting  of  copper  and  aluminum  foils  with  a  single 
mode  laser,  which  are  used  as  current  collectors  of  the  lithium- 
ion  battery  electrodes.  The  mathematical  model  includes  fluid 
flow,  heat  transfer,  evaporation,  multiple  internal  reflections,  free 
surface  evolution,  and  surface  forces  with  proper  material  phase 
changes.  Using  this  proposed  mathematical  model,  the  thresh¬ 
old  laser  parameters  for  cutting  as  well  as  penetration  time  and 
depth  have  been  investigated.  Given  the  presented  results,  we 
could  conclude  that  the  copper  laser  cutting  is  a  laser  intensity 
and  interaction  time  dependent  process.  The  aluminum  laser  cut¬ 
ting  depends  more  on  laser  intensity  than  interaction  time.  The 
proposed  mathematical  model  is  validated  experimentally.  Exper¬ 
imental  and  simulation  results  of  pure  current  collector  materials 
show  a  strong  agreement  for  the  high-speed  cutting  cases.  A  small 
discrepancy  observed  for  low-speed  cutting  can  be  improved  by 
modifying  the  constant  absorption  coefficient  to  be  a  temperature 
dependent  absorption  coefficient.  A  bubble-shape  recast  observed 
from  experimental  results  could  deteriorate  the  cut  surface  of 


electrodes  for  lithium-ion  batteries.  Simulation  results  of  pure  cop¬ 
per  are  in  good  agreement  with  experimental  results  of  pure  copper 
with  single  side-coated  graphite  when  a  laser  scanning  speed  is 
3000  mm  s-1.  Experimental  results  of  pure  aluminum  with  single 
side-coated  LiCo02  match  well  with  simulation  results  of  pure  alu¬ 
minum  when  a  laser  scanning  speed  is  5000mms_1.  Taking  into 
consideration  active  electrode  materials  in  the  proposed  math¬ 
ematical  model  can  improve  the  simulation  results  for  current 
collector  materials. 

This  study  could  be  utilized  as  a  guideline  to  obtain  good  qual¬ 
ity  of  a  cut  surface  so  that  clear  cut  edge  without  reattachments 
and  bubble-shape  recasts  can  prevent  the  catastrophic  failure  of 
the  entire  module.  Therefore,  the  desired  quality  of  a  cut  surface 
can  be  obtained  without  defects  by  controlling  the  laser  power  and 
scanning  speed.  Furthermore,  an  understanding  of  the  influences 
of  the  laser  parameters  on  the  laser  cutting  for  the  current  collector 
materials  would  be  used  for  further  investigation  of  the  electrode 
cutting  behaviors. 
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